Measuring fl with galaxy streaming velocities 



P. G. Ferrei^al'^ M. Davis^ H. A. Feldman''■^ A. H. Jaffe^ R. Juszkiewicz*''''®''^, 
^ Theory Group, CERN, CH-1211, Geneve 23, Switzerland 
CENTRA, Instituto Superior Tecnico, Lisboa 1096 Godex, Portugal 
University of California, Berkeley, CA94720, USA 
* Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve, 

Switzerland 

^Department of Physics and Astronomy, University of Kansas, Lawrence, KS 6604-5 
^ On leave from Copernicus Astronomical Center, 00-716 Warsaw, Poland 
''Joseph Henry Laboratories, Princeton University, Princeton NJ 08544 

The mean pairwise velocity of galaxies, Vi2 has traditionally been estimated from 
the redshift space galaxy correlation function. This method is notorious for being 
highly sensitive to the assumed model of the pairwise velocity dispersion. Here 
we propose an alternative method to estimate Vi2 directly from peculiar velocity 
samples, which contain redshift-independent distances as well as galaxy redshifts. 
This method can provide an estimate of Cl^'^as'^ for a range of ag where £7 is 
the cosmological density parameter, while trg is the standard normalization for 
the power spectrum of density fluctuations. We demonstrate how to measure this 
quantity from realistic catalogues and identify the main sources of bias and errors 



1 A model of vi2{r) 

In this presentation we report on the the possibihty of using the "mean ten- 
dency of well-separated galaxies to approach each other" ll^ to measure the 
cosmological density parameter, il. The statistic we consider is the mean rel- 
ative pairwise velocity of galaxies, Ui2. It was introduced in the context of the 
BBGKY theory El, describing the dynamical evolution of a collection of parti- 
cles interacting through gravity. In this discrete picture, Vi2 is defined as the 
mean value of the peculiar velocity difference of a particle pair at separation 
In the fluid limit, its analogue is the pair-density weighted relative velocity 



1 + U^) 

where and Sa = Pa/{p) — 1 are the peculiar velocity and fractional density 
contrast of matter at a point r^, r = |ri — raj, and ^(r) = (6162) is the two-point 
correlation function. The pair-weighted average, (• • •)p, differs from simple 
spatial averaging, (■■•), by the weighting factor p^pa {PfP2)~^, proportional to 
the number-density of particle pairs. In a recent letterB, one of us has shown 
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that an excellent approximation to V12 is given by 

Vi2{r) = -|i/r/|(r)[l+a|(r)] , (2) 

C(r) - (3/r3) r ^(x) dx ^ |(r) [ 1 + ^r) ] (3) 
Jo 

Here a is a parameter, which depends on the logarithmic slope of ^(r), while 
/ = d InD/dlna, with D(a) being the standard linear growing mode solution 
and a - the cosmological expansion factor. Finally, H = 100 h km s~^ Mpc 
is the present value of the Hubble constant. This approximate solution of the 
pair conservation equation, accurately reproduces, results of high resolution 
N-body simulations in the entire dynamical range B. 

If we restrict ourselves to r = 10h~^ Mpc , one can use the APM catalogue 
of galaxies S for an estimate of galaxy correlation function, £, — {r/r — O)"'^ ; 
the slope at the separation considered is 7 = 1.75 ± 0.1 (the errors we quote 
are conservative). One obtains then 

v.^ilOh-^ Mpc) = - 605 CT8^rj"-6 (1 + 0.430-82) /(I + 0.38cr82)2km/s. (4) 

The above relation shows that at r = 10h~^ Mpc, U12 is almost entirely deter- 
mined by the values of two parameters: and fl. It is only weakly dependent 
on 7. The uncertainties in the observed 7 lead to an error in Eq. ^of less than 
10% for CTs < 1. 



2 The estimator 

Since we observe only the line-of-sight component of the peculiar velocity, s 4 = 
Ta ■ Va/t = fji • Va, rather than the full three-dimensional velocity Va, it is not 
possible to compute V12 directly. Instead, we propose to use the mean difference 
between radial velocities of a pair of galaxies, ( Si — S2 )p — W12 f ■ (fi + ^2)72, 
where r = "Fi— r^. To estimate V12, we use the simplest least squares techniques, 
which minimizes the quantity x^(^) = J2a b H^a — Sb) — Pab Vi2{r)/2]'^ , 
where Pab = f ■ {^a + r^) and the sum is over all pairs at fixed separation 
r = |r_4 — r^l. The condition dx^/ dvi2 = implies 

I^Pab 

The above expression is a sum over positive quantities and so is stable. This 
estimator is appropriate to be applied to a point process which will sample an 



2 



underlying continuous distribution. The sampling is quantified in terms of tlie 
selection function, <p{r). The continuum limit of Eq. ^ is then 

~ j-^-j _ 2 / dm^ dm^ $12 (fii - 32)^12 
/ drrii dm2 $12^12^ 

with drriji — d^fA, and a two-point selection function given by $12 = 
i^ud'^i — ^2]" r) (t>{ri) (j){r2) , where Sd is the Dirac delta function. For ease of 
notation we shall denote the denominator in Eq. |^ by W{r). If we take the 
ensemble average of Eq. ^ we find that ('5i2('')) = Vi2{r)- Note that, unlike 
the estimators for the velocity correlation tensor proposed inci, the ensemble 
average of the estimator is Ui2(?') independent of the selection function. For 
an isotropic selection function this estimator is insensitive to systematic effects 
such as a bulk flow, large scale shear and small scale random velocities (as one 
might expect from virialized objects). 



3 Biases and errors 



How unbiased is this estimator? We have applied our statistic to mock cata- 
logues extracted from N-body simulations of a dust-filled universe with = 1 
and P{k) oc 1/fc. 

In Figure |l|(a) we plot 'Di2(f) with one standard deviation calculated with 
20 mock catalogues extracted with a deep selection function. Each catalogue 
has a different observation position within the simulation volume and so an 
average over this set should resemble a true ensemble average. The mean is 
consistent with what one would expect from a direct calculation with Eq. |^ 
(which is plotted in Figure |l|(a) as a solid line) . We have also performed this 
analysis without collapsing the cores; the results changed by very little. 

We repeat this calculation for a set of 9 catalogues all constructed from the 
same observation point using a deep (Figure |i|d) or shallow (Figure Qc) selection 
function to randomly sample a fraction of galaxies within the simulation box. 
The variance in 1)12(7') is now solely due to finite sampling ("shot noise"); for 
catalogues with 2000 to 3000 galaxies we expect the variance to be \/2~y/3 
times larger. We find that a shallow selection function changes the functional 
form, or slope, of the mean, making it a more sharply decreasing function of r 
than the ensemble average. It is therefore crucial when analyzing a catalogue 
to restrict oneself to scales much smaller than the effective cutoff scale of the 
selection function. 

Errors in distance measurements will naturally affect our results; the best 
estimators use empirical correlations between intrinsic properties of the galax- 
ies and luminosities and lead to log-normal errors in the estimated distance of 
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Figure 1: Vi2{r) (points) and its variance (dashed lines) evaluated from mock catalogues: A. 
random observers with a deep selection function compared to Vi2{r) evaluated from Eq, 
(solid line); b) A fixed observer with a deep selection function; c) a fixed observer with a 
shallow selection function. The variance is estimated from the scatter over 20 (a) or 9 (b,c) 

mock catalogues. 



around 20%. These errors will naturally lead to biases in cosmological estima- 
tors involving distance measurements and peculiar velocities and are generi- 
cally called Malmquist bias. We shall model our errors assuming a TuUy-Fisher 
law which resembles that inferred from the Mark III Ccitalogue. To correct for 
Malmquist bias we use the prescription put forward incl. 

In Figure H(a), we plot the results for the uncorrected simulations; Malmquist 
errors systematically lower the values of on small scales while enhancing its 
amplitude on large scales (where the effect should be more dominant). How- 
ever in Figure |2|(b) we show that with the correction for general Malmquist 
errors to the distance estimator, it is possible to overcome this discrepancy. 
The l-fj errors now encompass the true Vi2 over a wide range of scales. 

4 Conclusions 

In this contribution we report on a recent proposal to estimate the mean pair- 
wise streaming velocities of galaxies directly from peculiar velocity samples. 
We identified three possible sources of systematic errors in estimates of Vi2 
made directly from radial peculiar velocities of galaxies. We also found ways 
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Figure 2: Vi2{r) and its variance evaluated from 100 mock catalogues with errors (described 
in the text) and the full selection function. The solid points are the Vi2 of the error-free 
simulation seen from the same observation point, the solid line is the mean and dashed lines 
are the 1 cr. a) uncorrected distances; b) distances corrected for Malmquist bias 

of reducing these errors; these techniques were successfully tested with mock 
catalogues. The potential sources of errors and their proposed solutions can 
be summarized as follows. 

(1) On the theoretical front, assuming a linear theory model of Vizir) at 
r « lOh^^ Mpc can introduce a considerable systematic error in the resulting 
estimate of o'g'^il*''^. For example, if CTg = 1 using the linear prediction for Vi2 
at r = lOh^^ Mpc would introduce a 25% systematic ermr (see eq. Q). We 
solve this problem by using a nonlinear expression for Wiau 

(2) On the observational front, a shallow selection function induces a large 
covariance between on different scales. This must be taken into considera- 
tion by measuring {'12(7') only on sufficiently small scales. A rule of thumb is 
that for estimating at lOh^^ Mpc, the selection function should be reason- 
ably homogeneous out to at least 3Qh~^ Mpc. 

(3) Finally, care must be taken with generalized Malmquist bias due to log- 
normal distance errors; these induce a systematic error in Vi2- We have shown 
that, under certain assumptipns about selection and measurement errors, the 
method of Landy & SzalajO for corrected distance estimates allows one to 
recover the true Vi2. Naturally, this particular correction must be addressed on 
a case-by-case basis, given that different data sets will have different selection 
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criteria and correlations between galaxy position and measupanent errors. _ 
In a future publication we shall analyze the Mark III H and the SFI El 
catalogues of galaxies with this in mind. 
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